Using jet mass to discover vector quarks at the LHC 
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We illustrate the utility of jet mass distributions as probes of new physics at the LHC, focusing 
on a heavy vector-quark doublet that mixes with the top as a concrete example. For 1 TeV vector- 
quark masses, we find that signals with greater than 5<r significance can be achieved after 100 fb -1 . 
More generally, jet mass distributions have the potential to provide signals for heavy states that 
produce highly boosted weak gauge bosons and/or top quarks. 
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I. INTRODUCTION 

Various extensions of the Standard Model include ad- 
ditional quarks in vector representations of the standard 
model gauge group. Examples include little Higgs theo- 
ries [H, M [1] , topcolor models [4] , and extra dimensional 
theories with bulk fermions. Vector-like matter is also 
invoked in some supersymmetric models, for instance 
as the messengers of supersymmetry breaking in gauge- 
mediated models [j| . If light enough, vector-like quarks 
would be produced copiously at the LHC, and the details 
of how they might be discovered would depend on how, 
or whether, they decay. 

In this paper we concentrate on heavy quarks that de- 
cay into gauge bosons and top quarks. In particular, we 
will be interested in the case where the quarks are so 
heavy that the W's, Z's, and tops which they produce in 
turn yield highly collimatcd decay products that cannot 
typically be resolved into separate jets. The invariant 
masses of individual jets then become potentially useful 
quantities to study when attempting to pick out signals. 
Jet mass was originally used as a probe of QCD [6j , and 
has since been shown to be useful in studies of clastic 
WW scattering at high energy Q- Here we employ jet 
mass to provide signals for vector-quarks, and show that 
its usefulness persists after detector effects are taken into 
account. It is clear that jet mass could also be helpful in 
other collider studies with similar kinematic properties. 

One well-studied scenario with vector-quarks features 
electroweak singlets T + T, with hypercharge ±2/3. This 
is the extra fermion content in the littlest Higgs model 
If T mixes with the the top quark, then in the 
regime where the vector-quark mass is much larger than 
the top mass, the branching ratios for the decays of T 
are approximately Br(T -> bW + ) ~ 2Br(T -> tZ) ~ 
2Br(T — > th) ~ 1/2. Moreover, provided that the mix- 
ing is large enough, the T quarks can be produced singly 
by i-channel exchange of a W boson, with a b quark in 
the initial state 0] , the importance of the this being that 
the cross section for single production falls off less dra- 
matically with increasing T mass than the cross section 
for TT pair production. In Q, it was estimated that for 
order-one mixing, the discovery reach after an integrated 
luminosity of 30CT 1 fb at the LHC is ~ 2-2.5 TeV. On 
the other hand, if the mixing is small, or if the T quark is 
light enough, then QCD pair-production dominates. This 



case was studied in [l(J , where it was estimated that the 
discovery reach after 30CP 1 fb at the LHC is ~ 1.1 TeV. 

Here we study the case where the vector-quarks are 
instead electroweak doublets Q + Q, with hypercharge 
±2/3. Electroweak doublet vector-quarks appear, for in- 
stance, in the little Higgs model of and the topcolor 
model of [ll[ . We will argue below that it is reasonable to 
imagine that the decays of T and B, the upper and lower 
components of Q, are induced by the mixing of T with 
the top quark. In this case, B decays to tW~ , just as if 
it were a fourth-generation down quark. The prospects 
for discovering such a particle at the LHC have been ex- 
plored in [T3] • The approach taken there is to search for 
W candidates by looking at dijet invariant masses, find- 
ing top candidates by looking at the invariant masses of 
the the W candidates and ^-tagged jets, and finally, look- 
ing at the invariant mass distribution for the W and top 
candidates. The peak in this distribution ends up being 
rather broad, so that even for a B mass of 640 GeV and 
an integrated luminosity of 100 -1 fb, separation of signal 
from background looks challenging. 

Our strategy, which we will apply to TeV-mass -B's, 
will be to find signals in jet mass distributions. We will 
see that these distributions feature bumps around mw 
and m t , and that the invariant mass distribution for the 
candidate top and W jets is peaked near the vector-quark 
mass. 

In the next section we outline the model and list our 
assumptions about its free parameters. The main anal- 
ysis, requiring at least one high px lepton, is presented 
in Sec. III. There we list the relevant backgrounds and 
propose cuts that give a convincing signal. A dilepton 
analysis is presented in Sec. IV. 



II. THE MODEL 

Our vector quark doublet Q + Q has a mass term and 
also Yukawa couplings with the third generation quarks, 

£ maS s = MQQ + {X u Qu%h + X' u q 3 u%h + 

X d Qd c 3 h + X' d q 3 d c 3 h + h.c). (1) 

Electroweak symmetry breaking induces mixing of the 
heavy quarks with the third generation. We neglect mix- 
ings of the vector quark with the lighter two generations. 
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The main motivation for doing so is to simplify the anal- 
ysis, but the quark mass hierarchies themselves make this 
a reasonable starting point. Moreover, the mixings with 
the lighter generations are more tightly constrained ex- 
perimentally. Without loss of generality, we have defined 
Q to be the single linear combination of Y = 1/6 dou- 
blets that couples directly to Q, so there is no q% Q mass 
term. 

If we adopt the reasonable assumption that X u and A^ 
are of comparable size, and so are Xd and A^, then in the 
heavy vector-quark regime that interests us, M ^ m t , 
we have 

m t ~ X' u v, mt ~ X' d v, ~ tub — M, (2) 

where my and mg are the masses of T and B, the up 
and down components of Q. The ratio of the couplings 
of Q to the bottom and top quarks is then approximately 

(X d /X' d ) x (A'„/A„) x KM). (3) 

We will assume that this product is sufficiently small 
that decays of B and T directly into bottom quarks can 
be neglected. This assumption is motivated by the fact 
that mt, is much smaller than m t , and also by the fact 
that mixing with the bottom is constrained by Z-pole 
data. Mixing with the top is less constrained, as the 
electroweak couplings of the top quark have not been 
measured precisely. In this case, the equivalence princi- 
ple tells us that in the large M limit, the branching ratios 
of B and T are 

Br(B -> tW-) a 100% (4) 

Br(T -> tZ) ~ Br(T -> th) ~ 50% (5) 

Thus B decays just as if it were a fourth-generation down 
quark, while T has two possible final states, and all de- 
cays produce top quarks. 

Although electroweak symmetry lifts the degeneracy 
between B and T, the splitting S is only 

6 = m T -m B c 15 GeV x f^x (6) 

and the rate for the three-body decays T — > B ff is pro- 
portional to S 5 /v 4 . Even neglecting phase space factors, 
the ratio of this rate to the two-body rates is proportional 
to (rrit/M) 6 , and so we will neglect the three-body decay 
entirely in what follows. 

III. SINGLE LEPTON ANALYSIS 

We now explore the ability of the LHC to probe this 
model through final states with at least one lepton. For 
concreteness, we will fix M to be 1 TeV. Our discussion 
will focus on the production and decay of B particles, as 
our method will be far more sensitive to £Ts than to T's. 



For our calculations we will take the Higgs mass to be 
120 GeV, although our results are not very sensitive to 
its value. 

The heavy B can be produced singly in association 
with a top quark, with a cross section that depends on 
the heavy quark mass and on the amount of mixing. 
For large M, the cross section is roughly proportional to 
(\ u v) 2 /M 2 . Setting X u v = m t and M = 1 TeV, we find 
using Madgraph [13[ that the leading-order cross section 
for single B/B production is 14 fb. By comparison, the 
pair production cross section, which depends only on the 
heavy quark mass, is 60 fb at NLO with gluon resum- 
mation In what follows, we will focus on the pair 

production process exclusively. It is possible that single 
production may allow discovery of heavier i?'s if the mix- 
ing angle is large enough; on the other hand, if the mixing 
angle is much smaller than m t /M, single production is 
not likely to be of any help at all. 

In our analysis, we use MadGraph/MadEvent [l3[ to 
generate signal events at the parton level, taking the 
renormalization and factorization scales to be twice the 
heavy quark mass, and rescaling the cross section to agree 
with the NLO result. These are passed to Pythia 6.325 
fl5| to simulate initial and final-state radiation, multiple 
interactions, and hadronization. We use Alpgcn 2.06 fl6j ] 
and Pythia 6.325 to generate background events, and to 
obtain inclusive event samples we apply the MLM pre- 
scription for jet-parton matching [L 71 ] as implemented in 
Alpgen 2.11. The one exception is that jet-parton match- 
ing is not performed for single-top processes. For all cal- 
culations we use the PDF set CTEQ5L 

We pass showered events to the PGS-4 detector sim- 
ulator [l9[. Since we arc interested in jet mass, for our 
purposes the most crucial settings of the detector simu- 
lator are the energy resolution of the hadronic calorime- 
ter, which we take to be AE/E = 0.8/^E/GcV, and 
the granularity of the calorimeter, which we take to be 
A(j) x Ar] = 0.1 x 0.1. The PGS detector simulator uses 
the kr algorithm for jet clustering (2(j. We adopt refer- 
ence cone size AR cone = 0.5, and set the threshold trans- 
verse energy for a cell to be included in the clustering at 
5 GeV. This large threshold leads to a slight underesti- 
mation of jet energies, but turns out to be helpful for 
reducing the degradation of the signal from multiple in- 
teractions. We use the heavy-flavor tagging efficiencies 
included in the PGS code, which are based on the re- 
sults of a vertexing algorithm applied to CDF calibration 
data. Finally, the lepton isolation criteria are as follows: 
the total pt of tracks within AR = 0.4 of the lepton is 
required to be less than 5 GeV, and the total Et in the 
3x3 calorimeter grid with the lepton's cell at the center 
(excluding the Et of that central cell) is required to be 
less than 0.1 or 0.1125 times the Et in the central cell, 
for electrons and muons, respectively. 

To suppress backgrounds, we adopt the following cuts: 

• ^t+J^Pt > 1800 GeV, where the sum is over all 
photons, leptons, and jets having pt > 20 GeV 
and |rj| < 2.5 
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• at least one lepton (e. fi) with pr > 100 GeV and 
M < 2.5 

• j/ T > 100 GeV 

• at least one 6-tagged jet with pr > 20 

• Ai?y > 1.0, where Ai?y is the separation be- 
tween the hardest lepton and the closest jet having 
Pt > 20 GeV. This cut is useful for reducing the ti 
background, in which the leptons produced by the 
highly boosted tops are typically quite close to the 
b quarks. 

• St > 0.1, where St is the transverse sphericity. 
Given the 2x2 tensor SV, = ^aPfPj 1 ' where i 
and j label the two directions perpendicular to the 
beam, and a labels the jets, leptons, and photons 
having pt > 20 GeV and |r/| < 2.5, along with 
the missing pt, we define St as twice the smaller 
eigenvalue divided by the trace (so < St < 1 is 
always satisfied). 

After these cuts, the dominant backgrounds are W+ jets, 
ti, and Wbb, followed by WW and tW . We will assume 
that with these cuts we can neglect the background from 
QCD with a fake lepton, although to seriously address 
this background source would presumably require a study 
of fake rates using LHC data. 

To obtain large enough background samples, we im- 
pose the generator-level (pre- showering/hadronization) 
cut on all background processes except for single-top pro- 
duction: 

p't+Y.Pt > 1500 GeV, (7) 

where the sum is over the final-state partons (e.g. t, i, 
and any extra light quarks or gluons for the ti+ jets sam- 
ple) . We estimate how much of the background we lose 
by generating event samples below this cut, seeing what 
fraction pass the p't+^Zpt > 1800 GeV cut after de- 
tector simulation, and assuming that the rest of the cuts 
have roughly the same effect on the remaining events as 
for the sample that passed the generator-level cut. In this 
way we estimate that well under 10% of the background 
is neglected due to this generator-level cut. 

For the single-top final states tW, tq, tb, and tbW, 
multiple extra jets are not included at the matrix-element 
level, and we find that it is necessary to loosen this cut. 
After relaxing it to 

fc+^PT > 1000 GeV (8) 

for these processes, we again estimate that well under 
10% of the background is neglected. 

For all processes for which jet-parton matching is per- 
formed, we impose the following generator-level cuts on 
the extra light jets: pt > 20 GeV, \rj\ < 2.5, and a jet-jet 
separation ARjj > 0.7 (cuts of this nature are part of the 
jet-parton matching program). For Wbb, the minimum 



separation between extra light -jets and bottom quarks is 
also set to 0.7. The minimum cluster Et, rapidity range, 
and AR used for the jet-parton matching are then set to 
their Alpgen defaults of 25 GeV, \rj\ < 2.5, and 0.7, re- 
spectively. For tW jet-parton matching is not performed, 
and the extra jet in the iW+jet sample is required to have 
Pt > 20 GeV. 

For each background process, the factorization and 
renormalization scale Q is set to its Alpgen default, 
Q = ^2 m s m T i wnere the first term is the sum of the 
masses of any final-state gauge bosons, the second sum 
is over all final-state partons excluding the gauge bosons 
or their decay products, and where rrij, = p\ + m 2 , with 
m being the mass of the parton. 

In table |T] we list the numbers of events generated for 
the various backgrounds, the cross sections and corre- 
sponding integrated luminosities, and the numbers of 
events that pass all cuts. The cross sections are obtained 
by multiplying the Monte Carlo results by NLO K fac- 
tors. The factors we use are just rough estimates of the 
NLO effects, but neglecting them would certainly under- 
estimate the backgrounds. We take K — 1.5 for ti 21 1, 
K = 1.25 for W/Z + light jets ^%,K = 2 for Wbb 22 j, 
K = 1.1 for WW W, Z, and ZZ^LK = 1.1 for tq [H, 
K = 1.5 for tb 13, K = 1.2 for tW |25j|, and K = 1.5 for 
tbW, although we are not aware of NLO results for this 
last final state. In generating these samples, the states 
WW, tbW, and tW are decayed inclusively, the states it, 
ZZ, and WZ , are required to produce at least one lepton 
in their decays, and leptonic decays are required for the 
gauge boson or top quark in the W+jets, Z+jets, Wbb+ 
jets, tq, and tb samples. 

In table [II] we list the numbers of signal and back- 
ground events that pass the successive cuts, again scaled 
to 100 fb _1 . The lepton cuts are especially effective in 
reducing the ti background, because when the highly 
boosted tops decay leptonically, the lepton is often too 
close to the b quark to satisfy the isolation criteria. Jet 
mass distributions for the the signal and two of the most 
important background processes, W+jets and tt+jets, 
are shown in figure [T] Only jets with pt > 350 GeV are 
included, and distributions are normalized to 100 fb _1 
of integrated luminosity. The W + jets distribution de- 
creases significantly through the W mass, the ti + jets 
distribution is relatively flat in this region, while the sig- 
nal has a pronounced peak around the W mass. There 
is also a smaller bump around the top mass, although 
the ti + jets distribution also has a bump there due to 
the highly boosted tops. The presence of highly boosted 
W's also affects the ti+ jets distribution, but not nearly 
as dramatically as it does the signal. The distribution 
for Wbb+ jets, the other large background, steadily de- 
creases as the jet mass increases. 

A plot of the jet distributions for signal plus total back- 
ground and total background alone is shown in figure [2 
again including only jets with > 350 GeV. To estimate 
the significance of the peak around the W mass, we take 
three times the total number of jets in the 50-60 GeV bin 



4 





N 


cr(fb) 


£(fb" 1 ) 


N 


tt-\-3 jets (inc.) 


45,963 


516 


89.0 


80.8 


tt+2 jets 


11,333 


174 


65.3 


13.8 


tt-f-1 jets 


5,686 


86.4 


65.8 


1.5 


tt+O jets 


1,852 


34.0 


98.1 





W+4: jets (inc.) 


61,577 


725 


84.9 


102 


W+3 jets 


33,765 


375 


90.0 


27.8 


W+2 jets 


22,279 


237 


94.2 


19.1 


W+l jet 


9,081 


46.1 


197 





W+0 jets 


1,348 


5.37 


251 





Z+A jets (inc.) 


2,994 


75.8 


39.5 





^_l_3 j e t s 


6,126 


40.6 


151 


0.66 


Z+2 jets 


3,716 


27.0 


137 





Z+l jet 


2,550 


5.20 


490 





W/y[/_l_3 j e ts (inc.) 


9,471 


105 


90.1 


20.0 


WW +2 jets 


3,402 


42.4 


80.3 


8.7 


WW+1 jet 


1,145 


13.2 


86.5 





WW+O jets 


1,090 


3.08 


354 





ZZ+3 jets (inc.) 


212 


1.93 


110 





ZZ+2 jets 


77 


0.679 


113 


o 


ZZ+1 jet 


55 


0.262 


210 





ZZ+O jets 


125 


0.120 


1,041 





W(lv)Z(fJ)+3 jets (inc.) 


6,668 


26.3 


253 


9.1 


W(7*/)Z(f7 , )+2 iets 


3,703 


11.6 


321 


2.5 


W(MZ(/7)+l jet 


494 


3.50 


141 


2.1 


W(MZ(f7)+0 iets 


397 


0.421 


942 





VK(gg)Z(/ + r)+3 jets (inc.) 


1,355 


5.45 


249 





W(gg)Z(Z + r)+2 jets 


354 


2.24 


158 


0.63 


W(<jg)Z(Z+r)+l jet 


299 


0.762 


393 





W(gg)Z(Z + Z")+0 jets 


3,676 


8.12 x 10" 2 


4.5 x 10 4 





Wbb+2 jets (inc.) 


27,505 


177 


155 


102 


jet 


718 


12.0 


59.8 


13.4 


V(/b&_l_0 jets 


446 


1.53 


291 


1.4 


tW+1 jet 


20,000 


335 


59.7 


36.8 




17,771 


78.1 


228 


12.3 


tg 


5,487 


51.9 


106 





46 


950 


9.72 


97.7 







2,400 


60.6 


39.6 


15.2 


BB 


50,479 


60.0 


841 


210 


TT {HZ) 


7,951 


30.0 


265 


19.2 


TT (ZZ) 


7,954 


14.1 


564 


14.9 


TT (HH) 


7,969 


16.0 


498 


6.8 



TABLE I: For signal and background processes, the numbers of events generated N ge n, cross section, corresponding integrated 
luminosity, and number of events that pass all cuts N pas s, rescaled to an integrated luminosity of 100 fb _1 . For the background, 
the generator-level cuts described in the text are imposed. For the samples labeled "inc.", extra jets from showering are allowed 
when the MLM prescription for jet-parton matching is applied. The cross sections in this table include K factors to approximate 
NLO effects, as described in text. 
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signal(B) 


tt+ jets 


W+ jets 


Wbb+ jets 


generated 


6,000 


80,995 


138,801 


19,053 


J^Pt > !800 GeV 


2,610 


21,272 


44,175 


6,197 


lepton p T > 100 GeV 


864 


2,791 


12,634 


1,548 


j/ T > 100 GeV 


745 


2,035 


8,857 


1014 


at least one b-tag 


387 


1,009 


483 


302 


ARij > 1.0 


246 


182 


314 


210 


S T > 0.1 


210 


96 


149 


117 



TABLE II: For the B portion of the signal and the dominant 
background processes, the numbers of events that pass the 
successive cuts, scaled to an integrated luminosity of 100 fb~ . 
For the background processes the first row gives the number 
of events after the generator-level cut described in the text. 
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Jet Mass (GeV) 




100 150 200 
Jet Mass (GeV) 



300 



FIG. 1: Jet mass distributions for jets with pr > 350 GeV, 
for events that pass the cuts described in the text. We take 
100 fb _1 for the integrated luminosity. 



(187) as a background value to compare with the total 
number of jets in the 60-90 GeV bins (281), giving a 6.9a 
excess. More conservatively, taking the total number of 
jets in the 30-60 GeV bin (218) as the background value 
gives a 4.3c excess. Finally, taking the total number of 
jets in the 40-70 GeV bin (200) as the background value 
for the total number of jets in the 70-100 GeV bins (280) 
gives a 5.7 a excess. For each of these three measures, the 
standard model contribution to the number of events in 
the peak is smaller than the standard model contribution 
to the estimated background value. 

The PGS detector simulator does not include particle 
deflection by the magnetic field, but to get a rough idea 
of how sensitive our results are to this effect, we follow 
[26j and impose a shift in azimuthal angle for charged 
particles in the signal samples, 



L (0.45/pi 



(9) 



where the sign of the shift depends on the charge of the 
particle. We find that our results arc not dramatically 
affected by this shift. The significance estimates above 
change to 6.8cr, 4.2er, and 5.9a, respectively. 

The T quarks do contribute somewhat to the signal, 
because their decays can produce Z bosons, which are 
not resolved from Ws using this method. However, this 



FIG. 2: Jet mass distributions for the signal plus total back- 
ground and for total background alone, for events that pass 
the cuts described in the text. As before, only jets having 
Pt > 350 GeV are included for each qualifying event, and we 
take 100 fb _1 for the integrated luminosity. 



contribution is relatively small. Recalculating the signifi- 
cance in each of the three ways described previously, this 
time including only the B contribution to the signal, we 
find excesses of 6.3er, 3.5a, and 4.5er, respectively. 

We have seen that the jet mass distribution for the sig- 
nal is peaked around the W mass and less so around the 
top mass, due to the presence of highly boosted Ws and 
tops. Because the B quark decays as B — > tW~ (and 
the T quark decays as T — > tZ half of the time), one 
might hope to observe a peak in the invariant mass dis- 
tribution of pairs of jets whose masses are near rnw and 
mt, respectively. So, for each event passing our cuts, we 




200 400 600 800 10001200140016001800 2000 
Invariant Mass (GeV) 



FIG. 3: Invariant mass distribution for pairs of W and top 
candidates, after 100 fb _1 of integrated luminosity. 

identify as W candidates all jets with masses satisfying 
\mjet ~ m w\ < 20 GeV, and we identify as top candi- 
dates all jets with masses satisfying |m jet — m t \ < 30 
GeV. Then, for each event we pair up the W candidates 
with the top candidates in all possible ways, and calcu- 
late the invariant mass for each pairing. A histogram of 



6 



the resulting distribution is shown in figure [3l 

One simple point is that this procedure strongly en- 
hances the ratio of signal to background. Beyond that, we 
see a clear peak near the heavy quark mass of 1 TeV. This 
procedure tends to give a peak shifted somewhat below 
the actual mass. After 100 fb _1 of integrated luminosity, 
fewer than ~ 20 pairs are obtained with invariant mass in 
the 900-1000 GeV bin, but a more significant peak would 
be achieved with greater luminosity, or perhaps simply 
by optimizing cuts and adjusting the jet mass windows 
used to identify candidate W's and tops. Alternatively, 
using jet mass in tandem with a more refined method 
for identifying tops might enhance the signal. In our 
analysis there are fewer top candidates than there are W 
candidates, so to increase the significance one would first 
concentrate on enhancing the top signal. 



IV. DILEPTON ANALYSIS 

The ratio of signal to background can be improved 
by requiring two leptons (this also leaves essentially no 
hadronically decaying gauge bosons in the it and WW 
backgrounds). In this analysis, we impose the same cuts 
as before, except that now we require at least two leptons 
with px > 20 GeV and \r)\ < 2.5, the hardest with px > 
100 GeV. We will present results with and without the 
&-tag requirement. We do not consider the background 
source W+ jets with a fake second lepton. 

Without the 6-tag requirement, 133 signal events re- 
main (103 from B production and decay), versus 104 
background events. The dominant backgrounds are 
WW+jets (40 events), it + jets (28 events), Z+jets (16 
events), and WZ+jets (14 events). The resulting jet- 




100 150 200 
Jet Mass (GeV) 



300 



FIG. 4: Jet mass distributions for the signal plus total back- 
ground and for total background alone, for events passing the 
cuts in the dilepton analysis. Only jets having pr > 300 GeV 
are included for each qualifying event, and we take 100 fb _1 
for the integrated luminosity. 



number of jets in the 30-60 GeV bin (66) as a background 
value for the total number of jets in the 60-90 GeV bin 
(102), we find a 4.4cr excess. 

After a b-tag is required, 72 signal events remain (54 
from B production and decay), versus only 12 back- 
ground events. The jet mass distribution, shown in figure 
[5l has a peak with greater than 5<7 significance. 
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100 150 200 250 300 

Jet Mass (GeV) 



FIG. 5: Same as figure|4j but with a b-tag requirement. 



V. CONCLUSIONS 

Jet mass distributions will provide a useful probe of 
new physics at the LHC. In particular, we have seen that 
they give signals in excess of 5cr for vector-quark doublets 
of TeV mass after 100 fb _1 of integrated luminosity. 

It would be interesting to try to use jet mass to test 
other models. There may be parameter space in su- 
persymmetric models for which methods similar to the 
ones we have used would be effective, e.g. if heavy 
charginos are produced copiously and decay dominantly 
to W bosons and LSP neutralinos. Another possible ap- 
plication is to warped-space models [27] |. For example, 
in ref. [HI the detection of Kaluza-Klein bottom quarks, 
which decay in the same way as the B quarks considered 
here, was considered for masses in the ~ 500 GeV range. 
For heavier masses the methods outlined here would be 
useful. In rcf. [29], the detection of Kaluza-Klein gravi- 
tons through their decays into gauge bosons was studied, 
and jet mass distributions might be helpful there as well. 
Finally, techniques for dealing with highly-boosted tops 
in the context of warped models have been proposed in 
[30| . It is possible that jet mass considerations could also 
help for this purpose. 



Note added 



mass distributions are are shown in figure [H this time 
keeping only jets with px > 300 GeV. Taking the total 



A week after we submitted this article to the hep-ph 
arxiv, Ref. [3l[ appeared. We would like to draw readers' 
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attention to this work, as it also uses jet mass as a probe 
of new physics. 



The work of DTS was supported in part by NSF grant 
PHY-0555421. 



Acknowledgements 

The work of WS was supported in part by the De- 
partment of Energy under Grant DE-FG-02-92ER40704. 



[1] 



[2 

[3; 

[4] 
[5] 

[6 

[7] 
[8] 

[9 
[10 

[11 

[12 
[13 



N. Arkani-Hamed, A. G. C ohen and H. Georgi , Nucl. 
Phys. Lett. B 513, 232 (2001) [arXiv:hep-ph/0105239| ; 
D. E. Kaplan and M. Sc hmaltz, JHEP 0310, 039 (2003) 
[arXiv:hep-ph/0302049l; H. C. Cheng and I. Low, JHEP 
0309, 051 (2003) |arXiv:hep-ph/0308199| . For reviews 
see M. Schmaltz and D. Tucker-Smith, Ann. Rev. Nucl. 
Part. Sci. 55, 229 (2005) [arXiv:hep-ph/0502182| , and 
M. Perelstein, Prog. P art. Nucl. Phys. 58, 247 (2007) 



Phys. B 529, 424 (1998) [arXiv:hep-ph/9801375] 



arXiv:hep-ph/0512128 
N. Arkani-Hamed, A. G. Coh en, E. Katz and A. E. N et 
son, JHEP 0207, 034 (2002) [arXiv:hep-ph/0206021| . 
I. Low, W. Skiba and D. Smith , Phys. Rev. D 66, 072001 
(2002) [arXivihep ph/0207243| . 

C. T. Hill, Phys Lett. B 266, 419 (1991 ); Phys. Lett. B 

345, 483 (1995) [arXiv:hep-ph/9411426| . 

M. Dine, A. E. Nelson and Y. Shirman , Phys. Rev. D 

51, 1362 (1995) [arXiv:hep-ph /9408384 ; Phys. Rev. D 

53, 2658 (1996) |arXiv:hep-ph/9507378| ' 

L. Clavelli, Phys. Lett. B 85, 111 (1979), Phys. Rev. D 

26, 91 (1982); L. Clavelli and H. P. Nilles, Phys. Rev. D 

21, 1242 (1980); L. Clavelli and D. Wyler, Phys. Lett. B 

103, 383 (1981). 

J. M. Butterworth, B. E. Cox and J. R. Forshaw, Phys. 
Rev. D 65, 096014 (2002) [arXiv:hep-ph/0201098l . 
T. Han, H. E. Logan, B. McE lrath and L. T. Wang , 
Phys. Rev. D 67, 095004 (2003) 
M. Perelstein, M 



arXiv:hep-ph/0301040 
E. Peskin and A. Pierce, Phys. Rev. D 
69, 075002 (2004) [arXiv:hep-ph/0310039l . 
G. Azuelos et al, Eur . Phys. J. C 39S2, 13 (2005) 
[arXiv:hep-ph/0402037] . 

J. A. Aguilar-Saavedra, Phys. Lett. B 625, 
234 ( 2005) [Erratum -ibid. B 633, 792 (2006)] 



arXiv:hep-ph/0506187 



Rev. D 64, 035001 (2001) 



>A. B. Popovic, Phys . 
|arXiv:hep-ph/0101123] . 
Atlas Collaboration Physics TDR, ATLAS T DR14. 
http:/ /madgraph.phys.ucl.ac.be/index.html| F. Mal- 

JHEP 0302, 027 (2003) 
T. Stelzer and W. F. Long, 
ommun. 81, 357 (1994) 
H. Murayama, I. Watan- 



toni and T. Stelzer, 



arXiv:hep-ph/0208156 



[14] 



Compu t. Phys. C 
|arXiv:hep-ph/9401258] : 
abe and K. Hagiwara, KEK-91-11. 

R. Bonciani, S. Catani, M. L. Mangano and P. Nason, 



[15] |http://www.thep.lu.se/staff/torbjom/Pythia.html 



T. Sjostrand, L. Lonnblad, S. Mrenna and P. Skands, 



[16] 

[17] 
[18] 

[19] 

[20] 

[21] 

[22] 

[23] 

[24] 

[25] 
[26] 

[27] 

[28] 

[29] 

[30] 

[31] 



|arXiv:hep -ph/0308153 
M. L. Mangano, M. Moretti, F. 
tau and A. D. Polosa, JHEP 



arXiv:hep-ph/0206293 . 



Piccinini, R. Pit- 
0307, 001 (2003) 



See e.g. http://mlm.web.cern.ch/mlm/talks/lund-alpgen.pdf 
H. L. Lai et al. [CTEQ Collaboration] , Eur. Phys. J. C 
12, 375 (2000) arXiv:hep-ph/9903282 . 
J. Conwayet al., http://www.physics.ucdavis.edu/ 
conway /research / software / pgs /pgs4-general .htm . 
S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock and 
B. R. Webber, Phys. Lett. B 269, 432 (1991). 
M. L . Mangano, M. Moretti , F. Piccinini and M. Trec- 
cani, arXiv:hep-ph / 061 1 129 

J. Campbell, R. K. Ellis and D. L. Rainwater, Phys. Rev. 
D 68, 094021 (2003) [arXiv:hep-ph/0308195| ; C. Oleari 
and D. Zeppenfeld, Ph ys. Rev. D 69, 093004 (2004) 
[arXiv:hep-ph/0310156 



Oleari and D. Zeppen- 
B. Jager, C. Oleari and 
D 73, 113006 (2006) 



J- „ 

G. Bozz i, B. Jager, C. 
feld, |arXiv:hep-ph/0701105| 

D. Zeppenfeld, Phys. Rev 
[arXiv:hep-ph/0604200j . 
J. Campbell, R. K. Ellis and F. Tramontano, Phys. Rev. 
D 70, 094012 (2004) [arXiv:hep-ph/0408158 . 

N. Kidonakis, [arXiv:hep-ph/0701080| ~ 

E. Richter-Was., D. Froidevaux and L. Poggioli, ATL- 
PHYS-98-131. 

L. Randall and R. Sundmm, Phys. Rev. Lett. 83, 3370 
(1999) [arXiv:hep-ph/9905221| . 

C. Dennis, M. K. Unel, G. Servant and J. Tseng, 



arXiv:hep-ph/0701158 



K. Agashe, H. Davoudiasl, G. Perez and A. Soni, 
|arXiv:hep-ph/0701186| 

K. Agashe, A. Belyaev, T. Krupovnickas, G. Perez 
and J. Virzi, |arXiv:hep-ph/0612015| A. L. Fitz- 
patrick, J. Kaplan, L. Randall and L. T. Wang, 
|arXiv:hep-ph/0701150| B. Lillie, L. Randall and 
L. T. Wang |arXiv:hep-ph/0701166 
B. Holdom, |arXiv:hep-ph/0702037 



